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1.  Introduction 


Nearly  three  decades  ago  Wannier  proposed  that  the  electronic  Bloch  states  of  a  crys¬ 
talline  solid  become  sets  of  new  states,  later  called  Wannier-Stark  (WS)  ladders,  when  a 
constant  and  uniform  electric  field  is  applied.(l]  This  has  generated  an  intensive  research 
field  with  lasting  interest.  There  has  been  some  controversy  about  the  very  existence  of 
such  states, [2]  but  it  has  been  made  clear  that  the  notion  of  WS  ladders  is  valid  if  a  finite 
number  of  Bloch  bands  are  considered,  with  their  couplings  with  the  rest  of  the  states 
ignored.  [3]  The  criterion  for  ignoring  the  coupling  between  two  bands  is  that  they  are 
separated  by  a  gap  which  is  large,  such  that  the  process  of  Zener  tunneling  is  sufficiently 
slow  to  be  considered  unimportant  at  least  from  a  dynamical  point  of  view. [4]  Some  of  the 
basic  properties  of  the  WS  ladders  have  been  observed  experimentally.  [5] 

In  the  simple  case  where  a  given  Bloch  band  can  be  regarded  as  decoupled  from  the 
others,  the  nature  of  the  associated  WS  ladders  has  been  known  quite  well.  [6]  A  number 
of  discussions  have  been  made  regarding  the  choice  of  different  gauge  and  boundary  condi¬ 
tions,  all  revealing  the  basic  structure  of  the  frequency  spectrum;  uniformly  spaced  levels 
with  spacing  Fa/ft,  where  F  is  the  electric  force  and  a  the  lattice  constant  in  the  field 
direction.  In  a  static  representation  of  the  field,  the  eigenenergies  are  found  to  be  uni¬ 
formly  spaced  with  spacing  Fa,  and  the  eigenenergy  states  are  found  to  be  localized  more 
than  exponential  functions.  In  a  vector  potential  representation  of  the  field,  where  the 
Hamiltonian  is  no  longer  time  independent,  a  complete  set  of  solutions,  known  as  Hotiston 
functions,  is  found  for  the  time  dependent  Schrodinger  equation.  The  Houston  functions 
are  essentially  the  Bloch  functions  with  the  crystal  momentum  drifting  with  a  constant 
speed,  F,  in  the  field  direction.  The  Bloch  oscillation,  namely  the  periodic  recurrence  of 
the  the  Houston  fimctions  when  the  crystal  momentum  transverses  the  Brillouin  zone,  then 
gives  rise  to  the  uniformly  spaced  frequency  levels.  In  fact,  it  has  been  shown  that  any 
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initial  state  in  a  given  band  will  return  itself  in  a  period  of  Bloch  oscillation  (see  Luban  in 
ref  ,6). 

In  this  paper  we  study  the  effect  of  an  electric  field  on  an  isolated  Bloch  band  which  is 
split  by  a  potential  in  addition  to  that  of  the  crystal.  By  isolated  we  mean  that  the  spectral 
range  under  consideration  is  boimded  by  large  energy  gaps,  such  that  the  splitting  potential 
and  the  field  do  not  mix  it  significantly  with  the  rest  of  the  spectrum.  This  is  clearly  a 
situation  that  can  be  realized  in  an  experiment.  In  particulair,  we  think  of  a  super  lattice 
which  modulates  the  solid  in  one  of  the  crystal  directions,  with  the  electric  field  applied  in 
the  same  direction.[7]  As  far  as  the  motion  in  the  field  direction  is  concerned,  the  system 
can  be  ideally  modeled  by  an  one  dimensional  tight  binding  chain.  For  simplicity,  we 
will  take  the  hopping  amplitudes  to  be  nonzero  between  nearest  neighbors  only  and  to  be 
independent  of  positions.  The  splitting  potential  will  be  represented  by  a  set  of  position 
dependent  site  energies. 

We  wish  to  investigate  the  following  issues.  First,  we  would  like  to  show  analytically 
how  the  WS  ladders  of  the  parent  band  be  modified  under  a  weak  but  otherwise  arbitrary 
potential.  From  a  point  of  view  that  the  pieces  of  the  split  spectrum  are  going  to  be 
strongly  coupled  together  by  even  a  moderate  field,  the  solution  of  this  problem  should 
shed  some  light  on  a  strongly  coupled  multi-band  system.  Secondly,  when  the  potential  is 
periodic  (weak  or  strong)  but  not  constant,  the  parent  band  is  split  into  a  finite  number 
of  sub-bands.  In  this  case  it  is  known  that  there  will  be  n  sets  of  WS  ladders  under  a 
field,  where  n  is  the  number  of  sub-bands.  We  would  like  to  examine  how  the  energy 
(frequency)  levels  in  different  sets  are  placed  relative  to  one  another,  and  to  show  how 
their  positions  vary  as  functions  of  the  field  strength  or  other  parameters.  We  would  like 
to  see  the  behavior  of  the  ladders  in  the  extremes  of  weak  and  strong  fields,  and  to  see 
how  the  behavior  changes  from  one  extreme  to  the  other.  Also,  we  are  interested  to  know 
what  can  happen  when  a  pair  of  sub-bands  become  degenerate.  Thirdly,  about  the  time 
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evolution  of  the  system,  we  would  like  to  show  a  generalized  periodic  recurrence  of  the 
states,  which  has  been  expected  in  a  single  band  approximation.  Finally,  we  would  like  to 
develop  a  theoretical  framework,  which  exploits  the  symmetries  of  the  system  and  which 
allows  for  clean  and  easy  numerical  calculations. 

There  has  been  some  studies  of  multi-band  systems  in  the  literature,  but  many  of  them 
considered  the  case  of  weak  coupling  limit.[l,4,6,8,9]  It  is  in  this  limit  that  the  notion  of 
metastability  of  the  WS  ladders  has  been  extensively  investigated.  [9]  The  present  work  tries 
to  go  beyond  the  weak  coupling  limit,  and  to  emphasize  on  the  strong  coupling  regime. 
Fukuyama  et  al[10]  have  studied  a  tight  binding  model,  in  which  arbitrary  strength  of 
couplings  were  taken  into  account,  but  only  a  pair  of  bands  were  included.  There  have 
also  been  a  few  numerical  studies  on  finite  size  systems, [11]  but  a  clear  picture  has  not 
been  obtained  regarding  the  behavior  of  the  WS  ladders  as  functions  of  the  field  strength 
and  etc.  Finally  we  would  like  to  mention  a  branch  of  research  on  the  localization  of  the 
WS  states, [12]  where  the  regime  of  strong  coupling  has  also  been  dealt  with. 

Having  described  the  mun  scope  of  this  work,  we  now  proceed  to  outline  our  theoretical 
approach.  When  dealing  with  an  electric  field,  one  encounters  the  complication  that  it 
cannot  be  represented  by  gauge  potentials  which  are  uniform  both  in  spaure  and  time,  even 
if  the  field  itself  is  uniform.  Thus,  one  has  to  deal  with  a  Hamiltonian  which  is  either 
time  dependent  or  lacks  the  symmetry  of  spatial  translation  that  the  Hamiltonian  might 
have  in  the  absence  of  the  electric  field.  The  difficulty  has  been  reflected  in  the  earlier 
work  on  WS  ladders. [13]  It  is  well  known  that  the  WS  ladders  have  their  origin  in  the 
spatial  periodicity  of  the  crystal  potential,  so  it  would  be  nice  to  utilize  this  symmetry 
explicitly  in  one’s  treatment.  This  can  be  achieved  by  using  a  vector  potential  to  represent 
the  field,  but  the  Hamiltonian  then  becomes  time  dependent.  The  time  dependence  is  not 
a  problem  if  a  single  band  approximation  is  used,  but  it  can  cause  difficulty  in  a  more 
accurate  treatment  of  a  multi-band  system.  Earlier  utilization  of  the  spatial  symmetry  has 
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been  primarily  limited  to  inferring  some  general  structure  of  the  WS  ladders.  Therefore, 
when  calculating  the  spectrum  or  the  states,  one  often  took  the  static  representation  of  the 
field  and  used  fixed-end  boundary  conditions  or  alike.  The  boundary  effect  could  then  mix 
up  the  otherwise  perfect  structure  of  the  WS  ladders  and  make  it  difficult  to  determine 
how  different  sets  of  WS  ladders  place  themselves. 

This  embarrassing  situation  can  somehow  be  avoided,  if  it  is  realized  that,  when  mod¬ 
ified  by  appropriate  gauge  transformations,  the  space-time  translations  can  still  be  sym¬ 
metry  operations  of  the  Hamiltonian.  In  other  words,  one  has  to  replace  the  spau:e-time 
translations  by  their  gauge  covsiriant  forms  in  order  to  reveal  the  physical  symmetry  asso¬ 
ciated  with  the  uniformity  of  the  electric  field.  These  are  the  so  called  electric  trainslations 
originally  introduced  by  Ashby  and  Miller[14]  in  analogy  with  the  symmetry  of  magnetic 
translations.[15]  In  this  work,  we  will  utilize  the  symmetry  of  electric  translations  to  guide 
our  general  ideas.  The  actual  formulation  of  our  approach  will,  however,  be  made  in  an 
elementary  fashion. 

We  will  represent  the  electric  field  by  a  vector  potential  which  is  uniform  in  space 
and  linear  in  time.  We  will  take  a  periodic  boundary  condition  in  order  to  accommodate 
the  symmetry  of  spatial  translations  that  the  system  might  have.  The  Hamiltonian  is  no 
longer  time  independent,  so  we  will  focus  our  attention  on  the  evolution  operator.[16]  The 
symmetries  of  the  electric  translations  will  then  give  us  valuable  information  about  the 
frequency  spectrum  and  the  time  evolution  of  the  system. 

With  the  above  choice  of  gauge,  the  electric  translations  in  space  are  the  bare  trans¬ 
lations  in  space,  and  those  in  time  are  the  combined  time  and  momentum  shifts  (see  (2.4) 
and  (2.17)).  As  is  true  for  the  mi^netic  translations,  the  electric  translations  in  space  in 
general  do  not  commute  with  those  in  time,  unless  the  steps  of  spatial  and  time  translations 
satisfy  the  following  relation: 

6x6t  =  integer  x  h/F.  (1.1) 
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Due  to  this  restriction,  only  subgroups  of  the  electric  trainslations  can  be  diagonalized 
simultaneously  with  one  another  auid  with  the  evolution  operator.  The  subsequent  sections 
of  this  paper  are  basically  organized  according  to  three  classifications  of  the  symmetry 
subgroups,  each  emphasizing  a  different  aspect  of  the  system.  The  results  obtained  in 
different  ways  are  related  aoid  complementary. 

If  we  take  the  system  size,  Na,  as  the  basic  spatial  step  for  a  subgroup,  then  h/{FaN) 
is  the  finest  time  step  allowable.  The  consequence  of  the  symmetries  of  this  subgroup 
is  the  subject  of  study  in  section  2.  When  the  site  potential  is  periodic  with  period  pa, 
it  will  be  advantageous  to  consider  a  subgroup  with  pa  as  the  basic  spatial  step.  The 
finest  time  step  allowable  is  then  hl{Fap).  As  will  be  shown  in  section  3,  this  subgroup  of 
symmetry  operations  can  yield  a  great  deal  of  information  about  the  time  evolution  of  the 
system,  such  as  the  existence  of  a  set  of  invariant  subspaces  in  eixh  of  which  the  states  are 
localized  (more  than  exponential)  around  a  given  unit  cell.  The  WS  ladders  will  naturally 
manifest  themselves  through  the  periodic  structure  of  the  frequency  spectrum.  Section  4 
can  be  regarded  as  an  extension  of  section  3.  A  detailed  investigation  of  the  weak  and 
strong  field  limits  will  be  given  there.  The  weak  field  limit  is  in  fact  an  adiabatic  limit, 
where  one  can  see  how  the  Berry  phasesflT]  play  their  roles.  Some  numerical  results  for 
a  discrete  cosine  potential  will  be  given  to  illustrate  ideas  aind  to  cover  the  intermediate 
field  strengths.  In  section  5,  we  will  consider  yet  another  subgroup  of  electric  translations. 
The  basic  time  step  will  be  taken  as  hfFa,  which  is  the  period  of  Bloch  oscillations  of  the 
parent  band.  The  explicit  time  periodicity  of  the  Schrodinger  equation  will  be  exploited 
by  the  application  of  the  Floquet  theorem.  It  is  in  this  approach  that  weak  perturbations 
of  the  WS  ladders  of  the  parent  Bloch  band  can  be  conveniently  studied.  Examples  will 
be  given  for  a  single  impurity  potential  and  for  some  periodic  potentials. 

In  section  6,  We  will  summarize  our  results,  with  some  remarks  regarding  the  generality 
of  our  approach. 
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3.  Electric  translation  in  time  as  a  dynamical  symmetry 

In  the  simplest  form,  a  Bloch  band  is  described  by  the  tight-binding  Schrodinger  equa¬ 
tion 

Q 

0  =  -Rl2[rl>{n  +  1,0+  -  1, 01.  (2.1) 

with  a  periodic  boundary  condition 

il>(n  +  iV,  0  =  0,  (2.2) 

where  R  is  the  band  width,  and  N  the  total  number  of  sites.  A  constamt  phase  factor 
has  been  inserted  in  the  boundary  condition  for  generality  and  for  later  reference.  If  the 
system  is  perturbed  by  a  site  potential  Vn  and  an  electric  field,  equation  (2.1)  is  modified 
to  the  following  form 

*A^0(n,O  =  -iR/2(«*^‘*‘/^V'(»  +  1,0  -  1,0)  +  Vn^{n,t),  (2.3) 

where  F  is  the  electrical  force  and  a  the  lattice  constant.  The  electric  field  has  been 
represented  by  a  time-dependent  vector  potential.  This  introduces  a  time  dependence 
in  the  Hamiltonian,  and  therefore  renders  useless  the  convenient  characterization  of  the 
system  by  energy  eigenvalues  and  eigenstates.  It  is  known  that  the  time  dependence  can 
not  be  removed  by  any  choice  of  gauge  which  conforms  with  the  boundary  condition  (2.2). 
(Had  we  used  a  fixed-end  boundary  condition  (see  (2.22)),  the  time  dependence  cam  be 
removed  by  a  gauge  transformation  =  ejcp{iFant/h)f/f{n,t).  In  the  new  gauge,  the 

field  is  represented  by  a  static  potential  —Fan.)  The  motivation  of  insisting  on  a  periodic 
boimdary  condition  is,  of  course,  to  accommodate  the  symmetry  of  spatial  translation 
that  the  system  may  have.  It  will  be  seen  in  the  next  section  that  the  WS  ladders  can 
be  presented  most  cleanly  this  way.  Another  advantage  of  using  a  periodic  boundary 
condition  is  its  convenience  in  dealing  with  transport  problems.  We  would  like  to  explore 
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what  kind  of  dynamical  symmetries  still  remain  in  the  system,  assuming  at  the  moment  a 
general  static  potential. 

One  obvious  symmetry  is  the  time  periodicity  in  the  Hamiltonian,  which  is  associated 
with  the  Bloch  oscillation  of  the  parent  band.  A  number  of  useful  results  can  be  obtained 
by  utilizing  the  Floquet  theorem  applicable  to  lineu  differential  equations  with  periodic 
coefficients.  We  will  look  into  this  in  section  5. 

What  is  not  so  obvious  is  the  symmetry  of  the  following  type.  Suppose  ^(n,  t;  a)  is  a 
solution  of  (2.3)  and  (2.2),  then 

p/(n,  t;  a)  =  t  +  /r;  a),  /  =  integers  (2.4) 

are  also  solutions,  where  t  =  h/(FaN).  This  is  the  symmetry  of  a  subgroup  of  eleetrie 
translatioru  that  we  talked  about  in  the  introduction.  Because  of  the  boundary  condition 
(2.2),  the  basic  time  step  cannot  be  smaller  than  h/{FaN).  In  other  words,  the  boundary 
condition  implies  a  largest  possible  spatial  step  Na,  the  system  size,  such  that  the  sizes  of 
the  time  steps  have  to  be  bounded  from  zero  according  to  (1.1). 

It  is  convenient  to  introduce  the  evolution  operator  U{tft'’,a)  for  equation  (2.3)  with 
the  boundary  condition  (2.2).  By  definition,  if  ^(n,<)  is  an  arbitrary  solution  of  (2.3)  and 
(2.2),  then 

(2.5) 

n' 

The  symmetry  (2.4)  implies  that, 

91+1  («.  «)  =  <  +  »■;  a) 

=  ^  ^  a)g/(n,  t;  a). 

n' 

It  is  then  clear  that  a  state  ^(n,  t;  o)  in  the  first  time  interval,  0  <  t  <  r,  is  carried  over 
to  the  subsequent  intervals  by  repeated  application  of  the  operator  C/(t;a),  where 

o)  =  +  '■.‘1  “)■ 


(2.7) 


As  a  result,  the  system  is  completely  described  by  17(t,0;  a)  auid  U{t\  a)  for  0  <  t  <  r. 

We  now  proceed  with  a  general  analysis  of  the  operators  U  and  U.  Being  the  evolution 
operator  of  a  Schrodinr  equation  with  a  Hermitian  Hamiltonian,  U  is  unitary  and  satisfies 
the  following  relations: 

(2.8) 

where  t,  t'  and  t"  are  arbitrary  time  instants.  The  a  periodicity  of  the  boundary  condition 
(2.2)  implies  that 

U{tj'-,a-\-2nl{Na))  =  U{t,t';a).  (2.9) 

Moreover,  it  can  be  shown  that 


+  ha/F,t'  +  Aa/F;0)e-‘““", 


(2.10) 


which  is  true  for  any  Hamiltonian  that  is  independent  of  time  in  the  absence  of  the  electric 
field,  where  n  is  the  site  position  operator.  Using  these  general  relations  and  the  definition 
in  (2.7),  we  can  now  show  that 

=  +  r  +  fta/F,  t  +  htt/F;  0) 

=  +  r  +  ha/F,  r;  0)U(r,  0;  0)U(0,  t  +  ha/F;  0)  (2.11) 

=  U(t  -h  Aa/f’,0;2ir/(Wo))e‘2’""/^U(r,0;0)U(0,t  +  Aa/F;0) 

=  t/(t  +  Aa/F,0;0)U(0;0)U^{t  4-  ha/F,0;0) 

In  other  words,  0’(t;  a)  is  unitarily  similar  to  (7(0;  0),  and  therefore  has  a  set  of  eigenvalues 
independent  of  both  t  and  a. 

We  now  claim  that  the  characteristic  frequencies  of  the  system  are  of  the  form 

{u>j  +  2‘Knjl''y  Uj  —  integers,  j  =  1, 2, ... ,  AT},  (2.12) 
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where  {e  are  the  eigenvalues  of  0,  aind  it  is  understood  that  ufj  is  defined  modulo 

27r/T.  To  show  this,  we  make  an  expansion, 

ffo(n,fo;a)  =  52oj(to;oi)ttj(n,to;a),  (2.13) 

} 

where  Uj[n,tQ;  a)  is  the  jth  eigenstate  of  U{tQ;a).  Then  equation  (2.6)  yields 

9l{njQ]a)  =  ^aj{tQ\a)e~*^i^^Uj{n,tQ-,a)  (2.14) 

j 

In  terms  of  the  wave  function  V’l  this  is 

ili{n,tQ  +  Ir)  =  e“‘2’^”'/^^Oj(to;o)«~*"^’’^u;(n,to;ot)  (2.15) 

j 

The  overall  factor,  ^  gauge  factor,  and  will  be  canceled  out  in  the  expectation 

value  of  any  physical  observable.  It  is  easy  to  see  that  this  is  true  for  an  operator  that  is 
diagonal  in  n.  To  demonstrate  that  this  is  also  true  for  one  with  off-diagonals,  consider 
the  probability  current  from  site  n  to  site  n  +  1, 

+  1.0  -  +  1,1)1-  (2.16) 

When  t  is  increased  from  to  to  to  +  It,  the  gauge  factors  of  the  ^’s  in  the  above  expression 
are  cancelled  by  the  factors  produced  by  and  its  conjugate.  The  phase  feictors 

containing  the  cu’s  are  the  only  ones  left  which  depend  on  the  time  shift  Ir.  In  fact,  the 
physical  and  gauge  invariant  time-shifting  operator  should  be 

(2.17) 

instead  of  T(/r),  where  f'(Ir)  makes  a  pure  time  shift  of  It  on  a  state.  Therefore,  the 
possible  frequencies  of  the  system  are  indeed  of  the  form  in  (2.12). 


The  actir'  procedure  for  determining  the  eigen-frequencies  is  very  simple.  First,  we 
determine  0;0)  by 


{/(0,0;0)  =  1  . 
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where  the  boundary  condition  (2.2)  with  a  =  0  is  used.  The  matrix  is 

just  fy{t;oi)  with  zero  (  and  a.  The  eigen-frequencies,  being  independent  of  t  and  a,  are 
then  obtained  by  diagonalizing  U(0;0).  Here,  we  would  like  to  point  out  that  the  above 
procedure  may  be  carried  out  in  any  representation  which  is  convenient. 

The  upper  time  limit,  r,  of  the  time  dependent  problem  (2.15)  is  usually  small,  being 
inversely  proportional  to  the  system  size  N.  Therefore,  U  (r,  0)  cam  be  solved  with  ease 
and  high  precision.  When  N  is  very  large,  U (r,  0)  can  be  obtained  perturbatively.  To  first 
order  in  r,  we  have 

£/(r,0)  «  exp[-(*7h)  f  dtH{t)] 

JQ  (2.19) 

as  exp(— (ir/A)H(0)). 

In  this  limit  we  have  from  (2.7)  that 

£/(0;  0)  as  exp{  (0)  -  Fan}}.  (2.20) 

tfl 

Thus,  the  eigen-frequencies  of  O'  are  closely  related  to  the  eigenenergies  of  [^(0) —Fan],  the 
Hamiltonian  in  the  static  representation  of  the  field.  Our  approach  offers  the  convenience 
of  eliminating  the  boundary  effects  for  systems  of  finite  size. 

The  eigenstates  of  0^(t;a)  can  abo  be  easily  found.  We  first  look  for  the  eigenstates, 
{uj},  of  &(0;0).  Then,  according  to  the  last  equidity  of  equations  (2.11),  the  eigenstates 
of  i/(t;  a)  can  be  chosen  as 

+  Ha/F,  0;  0)  uj  .  (2.21) 

Before  we  close  thb  section  we  would  like  to  make  two  remaurks.  First,  the  dynamical 
symmetry  presented  above  should  be  applicable  to  more  general  systems  in  amy  dimension 
thain  our  one  dimensional  tight  binding  model.  Thb  fact  follows  from  the  generality  of  the 
symmetry  of  electric  tramslations. 

Secondly,  adl  the  results  will  remain  valid,  if  the  boundary  condition  (2.2)  is  replaced 
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by  a  fixed  end  type: 


V»(0,0  =V»(iV  +  l,0  =0  (2.22) 

Of  course,  the  group  of  dynamical  symmetry  in  this  case  can  be  extended  due  to  the  fact 
that  if  is  a  solution  of  (2.3)  and  (2.22),  then 

9/(n,t)  =  t  +  /r),  /  =  integers  (2.23) 

are  also  solutions,  with  r  arbitrary.  This  is  another  way  of  saying  that  the  time  dependence 
of  the  Hamiltonian  can  be  gauged  away  without  changing  the  boundary  condition  (2.23). 
Such  a  transformation  is  not  possible  with  the  periodic  boundary  condition  (2.2),  so  we 
had  to  make  the  restriction  r  =  hfFaN  in  (2.4). 

3.  Periodic  potential  and  WS  ladders 

In  this  section,  we  consider  the  case  of  a  periodic  potential  and  examine  the  conse¬ 
quences  of  the  symmetry  of  spatial  translation  on  the  frequency  spectrum  of  the  system. 

If  the  potential  in  equation  (2.3)  satisfies 

l^n+p  =  Vn  (3.1) 

and  the  system  size  iV  is  an  integral  multiple  of  the  period  p,  then  we  can  impose  the 
Bloch  condition  on  the  states  as 

i^(n  +  p,  t)  =  e‘*“PV'(n,  0  (3.2) 

where  0  <  k  <  2x/pa,  and  k  is  discretized  into  Nfp  possible  values  according  to  the 
boundary  condition  (2.2).  With  the  Bloch  condition  in  place  of  (2.2),  the  Schrodinger 
equation  (2.3)  needs  to  be  solved  only  for  n  in  a  period,  say  1  <  n  <  p.  The  arguments  and 
results  presented  in  the  last  section  can  then  be  transcribed  to  here  with  the  replacements: 
N  p  and  a  -*  k. 
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A 

Therefore,  if  f/(t,r;fc)  is  the  evolution  operator  of  equation  (2.3)  with  the  condition 
(3.2),  then  the  unitary  operator, 

U{t;k)  =  +  T,t-,k),  (3.3) 

will  have  p  eigenvalues  of  the  form  with  the  w’s  independent  of  t  and  k,  where 

T  =  h/Fap.  Moreover,  if  il/(n,t;k)  is  a  solution  of  (2.3)  and  (3.2),  then 

<  +  r;  *)  =  t).  (3.4) 

n* 

The  characteristic  frequencies  of  the  system  then  consist  of  p  sequences: 

{(jjj +2irnj{T,  nj  =  integers,  j  =  1,2,..., p},  (3.5) 

where  it  is  understood  that  the  cj’s  axe  defined  modulo  2)r/2'. 

The  average  frequency  spacing  in  the  spectrum  (3.5)  is  Fa/h,  the  same  as  in  (2.12), 
reflecting  a  conservation  of  level  density.  The  spatial  periodicity  of  the  system  has  imposed 
a  structure  in  the  spectrum  (2.12),  that  is,  a  periodicity  over  every  group  of  p  levels.  In 
other  words,  the  spectriim  has  been  organized  into  p  sets  of  levels,  with  a  uniform  spacing, 
Fapfh,  within  each  set.  These  correspond  exactly  to  the  WS  ladders  found  before  in 
various  ways.  In  this  sense,  we  say  that  the  system  is  characterized  by  p  sets  of  WS 
ladders,  each  represented  by  a  u;j. 

To  get  some  insights  into  the  behavior  of  the  states,  we  now  go  to  a  representation  with 
the  eigenstates  of  as  a  basis.  Consider  the  following  expansion  of  an  arbitrary  state: 

V'(«.  0  =  Cj{t',  k)uj{n,  t\  k),  (3.6) 

*  j 

where  uj  (n,  t;  A;)  is  an  eigenstate  of  0 (t;  A:)  and  satisfies  the  Bloch  condition  (3.2) .  The  wave 
function  inside  the  first  summation  symbol  in  the  above  expression  should  be  a  solution 
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of  (2.3)  and  (3.2),  and  therefore  is  connected  to  its  time  translation  as  in  (3.4).  Then  it 
follows  that 

0(n,t  +  IT)  =  ^  Jb)e-^i^^uy(n,t;  fc).  (3.7) 

i.* 

We  will  now  see  how  the  time  evolution  of  the  state  is  influenced  by  the  fact  that  the  u/’s 
are  independent  of  k.  It  will  be  convenient  to  consider,  for  a  sequence  of  time  (to  +iT,  /  = 
integers},  the  density  matrix  defined  by 

p(n,n';to  +  iT)  =  e‘2»"'/^V'(n,to  +  (3.8) 

where  the  gauge  invariant  form  of  the  time  shifting  operation  has  been  used.  (See  equation 
(2.17)  and  the  arguments  given  just  above  and  below  it  for  an  explanation  of  this  usage). 
Then  it  can  be  seen  that  the  partial  trace  of  the  density  matrix  over  jf, 

5]  fjjik,  (0  +  IT)  =  53  Cj(loi  t)c;((o;  *')  (3.9) 

J  3 

is  independent  of  the  time  shifting,  /T,  where  by  definition 

Pjj>(k,k';tn  +  IT)  =  ^5^uJ(n,to;fc)p(n,n';to  +  ^r)ttjv(n',to;*')>  (3-10) 

nn' 

with  {uj{n,tQ;k)}  normalized  over  an  unit  cell.  This  means  that  the  expectation  value  of 
any  physical  operator  acting  on  k  only  is  periodic  in  time.  One  trivial  example  of  such 
an  operator  is  the  crystal  momentum,  k,  itself.  A  more  interesting  one  is  the  projection 
operator  into  a  set  of  Wannier  functioiu  associated  with  the  Bloch  functions,  {ttj(n,  to\k)}, 
at  the  initial  time  More  specifically,  the  Wannier  functions  here  axe  defined  as 

wj{n  -  n^;to)  =  to;  *),  (3.11) 

where  labels  the  position  of  the  unit  cells  and  can  be  taken  for  its  value  as  0,  ±p,  ±2p, ..., 
etc.  The  arguments  following  (3.10)  then  indicate  that  the  subspace  spanned  by  {wj(n  — 
rifiito)]  j  =  1,  ..,p}  for  a  given  is  an  invariant  space  of  the  system  under  time  translations 
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by  IT.  In  the  ewe  of  p  =  1,  we  can  choose  as  the  eigenstate  of  U{tQ;k),  which  implies 
that  the  Wannier  function  will  be  strictly  localized  on  a  site.  Therefore,  a  state  initially 
localized  on  a  site  will  agaiin  be  localized  on  that  site  after  a  period  of  time  T  =  h/[Fa). 

We  now  examine  the  localization  property  of  the  Wannier  functions.  The  results  of 
last  section  (especially  (2.11)  and  (2.21))  indicate  that,  if  Uj(n,0;0)  is  the  jth  eigenstate 
of  &(0',0),  then  we  can  choose 

Uj{n,t-k)  =  ^ +  hfc/F,0;0)  Uj(n',0;0)  (3.12) 

n* 

as  the  jth  eigenstate  of  C{tQ;k),  where  the  summation  is  restricted  on  a  unit  cell.  The  n- 
independent  phase  factor,  ensures  the  periodicity  as  Jb  is  increased  by  a  reciprocal 

lattice  number,  27r/ (pa),  as  can  be  verified  directly  by  using  the  general  properties  of  U 
and  U  listed  in  (2.7)  through  (2.11)  in  the  last  section.  From  the  theory  of  initial  value 
problems  of  differential  equations,  (fnn'iH  +  /ljfc/F,0;0)  is  an  entire  function  of  its  time 
variable,  and  therefore  of  k.  To  summarize,  Uj{n,t]k)  can  be  chosen  such  that  it  is  a 
periodic  and  entire  function  of  A;,  and  that  it  satisfies  the  Bloch  condition  (3.2).  All  these 
properties  ensure  that  the  Wannier  functions  in  (3.11)  be  localized,  and  that  u;j  (n  — n^;to) 
go  to  zero  faster  than  any  exponential  function  as  |n  —  gets  large.  [19]  The  result  of  the 
last  paragraph  then  tells  us  that,  if  a  state  is  initially  localized  then  it  remaiins  so  in  time. 

4.  Weak  and  strong  field  limits  and  model  studies 

In  the  weak  field  limit,  the  Hamiltonian  in  (2.3)  varies  slowly  in  time,  and  it  can  be 
useful  to  go  to  a  representation  with  the  adiabatic  eigenstates  of  the  Hamiltonian  as  a 
basis.  Consider  the  expansion  of  an  arbitrary  state: 

(4.1) 

0,k 
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where  /?  is  a  band  label,  and  4>^{n,  a),  with  s  =  k  +  Ft/h,  is  an  adiabatic  eigenstate  of  the 
following  problem 

«)  =  — i2/2[e**“^^(n  +  1,  s)  +  e  —  1,  a)]  +  V’n^^(n,  a) 

(4.2) 

4>0{n  +  p,s)  =  <f>0{n,s) 

It  is  seen  from  the  above  equation  that  the  wave  function, 

«'*‘>"^^(n,a),  (4.3) 

is  just  a  Bloch  function,  with  wave  number  a  and  energy  c^(a),  of  the  Hamiltonian  in  (2.3) 
in  the  absence  of  the  electric  field.  It  is  well  known  that,  if  the  Bloch  bamds  do  not  cross 
one  another,  then  the  Bloch  function  can  be  chosen  as  a  smooth  and  periodic  function  of 
the  wave  number.  [19]  Therefore,  we  can  impose  the  aidditional  condition  on  ^  as 

^^(n,a  4-  2jr/pa)  =  e~’^’'"/P<^^(n,a).  (4.4) 

Now,  the  following  wave  function, 

Cpit;  fc)e‘^0^(n,  k  +  Ft/h),  (4.5) 

0 

in  the  expansion  (4.1)  should  satisfy  (2.3)  and  (3.2),  so  it  is  connected  to  itself  a  period  T 
time  later  by  C^(t;k)  as  in  (3.4).  It  then  follows  that 

C0(t  +  T;k)  =  J2^00'(t;k)C0,{t;k),  (4.6) 

0' 

where 

*)  =  E  *  +  Ftlh)tg.(n’,k  +  «/»)£?„„, (I;  k),  (4.7) 

nn' 

assuming  that  the  ^’s  are  normalized  in  the  unit  cell.  Thus,  U  is  unitarily  similar  to  U 
and  must  have  the  same  set  of  eigenvalues. 
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The  equation  of  motion  for  the  expansion  coefficients,  C^,  can  be  found  by  substituting 
(4.1)  and  (4.2)  into  (2.3),  with  the  results 

*)  =  €/j(*  +  Ft/K)Cg(f,  t)  +  5]  X;3I3,  (*  +  Ft/h)Cg,(,f,  k),  (4.8) 

0' 

where  the  inter-band  coupling  matrix,  X(s),  is  given  by 

X00'{s)  =  F^<f>}{n,a)i—il>i3,{n,s)  (4.9) 

n 

So  far,  the  results  have  been  exact  so  long  as  the  bands  do  not  cross  one  another.  The 
adiabatic  approximation  amounts  to  dropping  the  off-diagonal  elements  of  X(s).  Then, 
equation  (4.8)  can  be  easily  solved  to  yield 

*  T 

c^it  +  T;  *)  =  C^(t;  k)exp{^  (e^(*  +  Ft/h)  +  X^pik  +  Ft/h)]},  (4.10) 

with  the  evolution  operator  U[t]k)  diagonal  in  the  band  index,  and  with  the  eigen- 
frequencies  given  by 

=<  >  +  <  Xffff  >,  (4.11) 

where  the  angular  brackets  denote  aver^es  over  a  period  of  time,  T.  That  the  eigen- 
frequencies  are  independent  of  t  and  k  can  be  seen  from  the  fact  that  both  and  X^^i 
depend  on  t  and  k  through  the  combined  variable  [k  +  Ft/h)  and  are  periodic  in  time  with 
period  T. 

The  formula  (4.11)  says  that  hcjp  is  equal  to  the  average  band  energy  plus  an  extra 
term  <  X^p  >.  The  extra  term  is  in  general  nonzero,  and  can  be  related  to  the  Berry 
phase  associated  with  the  adiabatic  eigenvalue  problem  (4.2).  If  we  had  imposed  on  ^  the 
condition  of  parallel  transport, 

5^^J(n,s)^^^(n,3)  =0,  (4.12) 
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then  there  would  have  been  an  extra  phase  factor  multiplying  the  right  hand  side  of 
(4.4).  The  phase  is  caJled  the  Berry  phase  for  the  /9’s  adiabatic  level  of  (4.2),  with  its 
value  given  by 

^13  =  >.  (-'•13) 

where  <  >  should  be  evaluated  with  the  original  condition  (4.4). [17]  Since  the  Berry 

phase  is  independent  of  the  field,  the  term  <  Xpp  >  is  proportional  to  the  field.  It  is  seen 
that  (4.11)  is  really  the  first  two  terms  in  an  expansion  in  powers  of  the  field  strength. 

The  correction  to  the  adiabatic  approximation  (4.10)  will  be  exponentially  small  when 
~  */(?')»  is  small.  Since  Xppi  is  of  order  Fap/{2ir),  the  criterion  for  a 

good  approximation  is 

Fapl{2it^^^,)  C  1,  (4.14) 

where  is  the  minimum  gap  between  the  ^th  and  fi'th  energy  bands. [4]  The  corrections 
will  in  general  show  up  in  the  off-diagonal  elements  of  and  lift  the  degeneracies  (if  any) 
in  the  adiabatic  eigenvalues  This  fact  has  also  been  noticed  by  other  authors.(9] 

Now  consider  the  strong  field  limit.  We  can  then  make  a  first  order  time-dependent 
approximation  for  the  unitary  matrix  ^(7',0;0),  because  the  time  period  T  =  h/{Fap) 
will  be  small.  Therefore 

=  exp[(-iT/h)B^|, 

where 

^0(n)  =  sin(5r/p)[e‘’^/P^(n  +  1)  +  -  1)]  +  Vn^{n).  (4.16) 

2  TT 

The  corresponding  approximation  for  t^(0;0)  is  then 

i?„„-(0;0)  «  .-‘Wp  ^  (4.17) 

m=l 
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where  {€m»  Me  the  eigenenergies  and  eigenstates  (orthonormalized  over  a  unit  cell) 

of  (4.16).  When  T  -*  0,  the  matrix  (4.17)  will  be  diagonal  in  n,  due  to  the  completeness 
of  the  eigenstates  {^m(n)}.  The  p  diagonal  elements  are  all  different,  so  when 

the  off-diagonals  are  small  compared  with  the  minimum  level  distance,  2sin(»r/p),  the 
eigenvalues  of  U(Q;0)  can  be  approximated  by 

exp[-t(2irn/p  +  r?(n)/h)]  (4.18) 

where 

=  X)  «m^m(n)^;;,(»)  (4.19) 

m 

is  the  average  energy  on  the  site  n.  In  view  of  the  above  arguments,  the  condition  for  the 
strong  field  approximation  should  be 

«  1.  (4.20) 

We  will  now  study  a  concrete  model  to  illustrate  the  ideas  presented  above  and  in  the 
previous  sections.  Specifically,  we  take  the  site  potential  in  the  form 

Vn  =  V  cos(^n),  (4.21) 

where  q  and  p  are  mutually  inconmmensurate  integers  so  that  p  is  the  fundamental  pe¬ 
riod  of  the  potential.  Without  the  electric  field,  the  Schrodinger  equation  (2.3)  with  this 
potential  is  the  time  dependent  form  of  the  well  known  Harper’s  equation  or  the  discrete 
Mathieu’s  equation.  [20]  Harper’s  equation  has  been  extensively  studied  in  the  context  of 
two  dimensional  periodic  systems  in  a  uniform  magnetic  field.  Our  study  here  will  be 
applicable  to  such  systems  with  an  additional  electric  field  applied  on  them.[2l|  What 
we  like  to  know  particularly  is  how  the  eigen-frequencies  behave  when  the  field  or  other 
parameters  are  varied. 
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The  case  of  p  =  1  is  very  simple.  The  solution  of  (2.3)  and  (3.2)  can  be  found  exactly 
as 

V'fc(»»,t)  =  {-Rcc»(Fat'  jh  +  ka)  +  V”)].  (4.22) 

This  is  the  so  called  Houston  function  for  the  Bloch  band  under  consideration.  The  unitauy 
operator  U  is  now  a  scalar  and  can  be  read  off  from  the  above  as 

1  fhUFa) 

U  =  exp(—  dt  {-R coa{Fat'/h  +  ka)  +  V)] 

_  ^-i2irV/Fa 

—  c  , 

with  w  =  V/h,  which  is  independent  of  A;  as  expected.  The  single  set  of  WS  ladders  is  then 

{w  +  nFa/h,  n  =  integers}.  (4.24) 

We  expect  that  these  uniformly  spaced  levels  will,  in  generad,  be  modulated  to  non-uniform 
ones  in  the  presence  of  a  varying  potential;  the  modulation  will  be  periodic  if  the  potential 
is  periodic. 

When  p  >  1,  exact  solutions  are  hard  to  find,  and  numerical  calculation  has  to  be 
invoked.  We  have  considered  several  cases,  each  to  reveal  one  or  more  aspects  of  the 
behavior  of  the  eigen-frequencies.  In  the  figures  (larle)  presented  and  analyzed  below,  the 
reduced  eigen-frequencies  {fi(jjj/{Fa)j  j  =  l,...,p}  are  plotted  against  Rf{Fap),  namely 
the  inverse  field  scaled  by  the  hopping  amplitude,  R.  Since  wjT  =  ufjh/{Fap)  is  only 
defined  modulo  2)r,  we  found  it  convenient  to  fold  hufjl{Fa)  inside  the  range  [0.5,  p  + 
0.5).  The  reader  should  be  aware,  however,  that  the  whole  frequency  spectrum  extends 
periodically  both  upwards  and  downwards.  The  other  dimensionless  parameter  besides 
R/{Fap)  iaV/R.  We  have  taken  it  as  0.5  for  the  first  three  figures  and  0.1  for  the  last  two 
in  order  to  study  the  limit  of  small  band  splittings. 

Figure  la  corresponds  to  q/p  =  1/2.  In  the  strong  field  limit,  the  spacing  between  the 
two  b;’s  approaches  Fafti,  meaning  that  the  whole  spectrum  reduces  to  the  single  set  of 
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WS  ladders  (4.24)  of  the  parent  band.  The  initial  slopes  of  the  two  levels  agree  perfectly 

well  with  the  prediction  of  (4.18),  which  now  has  the  form: 

huJili^Fa)  =  1  -  2V/{Fap) 


(4.25) 


fUj2/{Fa)  =  2  +  2Vf(Fap). 

Deviation  from  the  above  behavior  becomes  pronounced  when  R/[Fap)  >  0.2.  As  the  field 
gets  weaker,  the  levels  approach  straight  lines  with  decreasingly  small  repulsions  at  their 
intersections.  The  straight  line  behavior  is  a  manifestation  of  the  adiabatic  limit.  (The 
band  structure  of  (4.2)  is  shown  in  figure  2a.)  The  expression  (4.11)  can  be  evaluated  with 
little  effort  to  give 


h(Jif{Fa)  =  g{RfFap) 


huj2l{Fa)  =  1  -  g{R/Fap), 


(4.26) 


where 


g  =  2  j  ds  yj (VIR)^  +  coe2(xs). 


(4.27) 


with  g  =  1.67761  d^t  V/R  =  0.5.  When  the  expressions  in  (4.26)  are  folded  in  the  range 
of  the  plot,  they  match  very  well  with  the  levels  in  the  figure  for  Rf{Fap)  >  1.  The  level 
repulsions  at  their  intersections  are  due  to  the  coupling  of  the  two  adiabatic  bands.  The 
size  of  the  gaps  generated  by  level  repulsions  decreases  exponentially  with  decreasing  field, 
that  is. 


^(/iu;/(J'o))  ~  exp(— AJ?/(Fop)),  (4.28) 

where  the  exponent  A  is  about  0.77  for  VfR  =  0.5. 

Figure  lb  is  for  q/p  =  1/4.  The  new  feature  appeared  in  this  case  is  that  two  of 
the  levels  are  horizontal  lines.  Explicit  calculation  shows  that  the  four  adiabatic  bamds  are 
symmetrically  located  about  e  =  0,  and  that  the  central  two  are  degenerate  (see  figure  2b). 
If  we  treat  the  degenerate  bands  as  one  big  band  and  apply  the  adiabatic  approximation 
to  it,  then  we  obtain  the  two  horizontal  levels  in  the  figure.  They  are  horizontal,  because 
the  average  energy  of  the  combined  band  is  zero.  The  Berry  phase  term  never  gives  rise  to 
a  nonzero  slope  in  our  plot,  and  it  only  shifts  the  levels  vertically  to  their  right  positions. 
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Figxire  Ic  corresponds  to  the  case  of  q/p  =  2/5.  As  noted  before,  the  levels  at  weaker 
fields  approach  those  of  the  adiabatic  approximation,  with  steeper  lines  corresponding  to 
bands  farther  away  from  the  zero  energy.  (See  figure  2c  for  the  band  structure.)  The 
adiabatic  level  corresponding  to  the  central  band  should  be  nearly  a  horizontal  line  at 
about  2.5  vertically,  but  is  interrupted  into  segments  due  to  the  repulsion  with  the  other 
adiabatic  levels.  On  the  other  hand,  the  initial  positions  and  slopes  of  the  levels  in  the 
strong  field  limit  sire  well  accounted  for  by  the  approximation  (4.18). 

In  the  zero  potential  limit,  the  levels  should  become  those  of  the  WS  ladders  of  the 
parent  Bloch  band,  that  is,  they  should  become  horizontal  lines  with  integral  heights  in 
our  plot.  To  study  the  effect  of  a  weak  periodic  potential,  let  us  consider  again  the  cases 
qjp  =  1/2  and  2/5,  but  with  V fR  =  0.1.  The  results  are  displayed  in  figures  Id  and  le. 
It  is  seen  that  the  levels  are  modulated  up  and  down  sinusoidally  about  the  WS  ladders 
of  the  parent  band,  as  the  field  strength  b  varied.  The  amplitude  of  modulation  increases 
with  decreasing  field,  while  the  frequency  of  modulation  approaches  to  a  constant.  A 
quantitative  analysis  of  these  results  will  be  given  in  the  next  section. 


5.  Floquet  theorem  and  weak  potential  limit 

In  this  section  we  will  make  use  of  the  symmetry  associated  with  the  explicit  time- 
periodicity  of  the  Schrodinger  equation  (2.3)  by  applying  the  Floquet  theorem.  [22]  As  one 
will  see  below,  this  approach  will  lead  to  a  very  compact  theory  for  the  WS  ladders  of  the 
parent  Bloch  band  under  weak  perturbation. 

The  theorem  tells  us  that  the  evolution  operator  of  equation  (2.3)  can  be  written  in 
the  form 
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(5.1) 


where  IV  vs  Hermitian  and  independent  of  time,  Md  Z[tyt')  is  unitary  and  periodic  in  t 
and  t'  with  period  /i/{Fo).  The  time  evolution  of  the  system  is  then  characterized  by  a  set 
of  frequencies  given  by  the  eigenvalues  of  W  plus  integral  multiples  of  Fafh.  Because  of 
its  additive  multiplicity,  the  frequency  spectrum  can  also  be  determined  by  the  eigenvalues 
of  the  following  operator 

^-iWh/iFa)  ^  u(^hl{Fa),0)  =  U,  (5.2) 

where  the  first  equality  comes  from  the  periodicity  of  Z  and  from  the  fact  that  Z  (0, 0)  is 
the  umty  operator.  By  definition,  U  is  just  the  evolution  operator  over  the  first  period. 
The  above  relation  indicates  that,  if  is  the  set  of  eigenvalues  of  U,  then  the 

frequency  spectrum  is  given  by  {u;m}  plus  integral  multiples  of  Fa/h. 

The  operator  U  here  also  has  a  simple  relationship  with  those  introduced  in  the  previous 
sections  in  the  same  notation.  It  is,  in  fact,  equal  to  the  Nth  power  of  that  in  section  2, 
and  the  pth  power  of  that  in  sections  3  and  4.  As  a  result,  the  frequency  spectrum  given 
in  this  section  will  correspond  to  those  modulo  Fa/h  in  the  previous  sections.  It  will  be 
helpful  to  notice  that  we  have  been  looking  at  the  problem  (2.3)  in  three  different  time 
scales:  h/{FaN)  in  section  2,  h/Fap  in  sections  3  and  4,  and  h/Fa  here  in  this  section. 
The  degeneracy  of  the  frequency  spectrum  when  viewed  from  a  longer  time  scale  should 
correspond  to  a  periodic  structure  of  the  spectrum  viewed  from  a  shorter  time  scale,  and 
vice  versa. 

When  the  potential  Vn  is  weak  and  the  field  is  strong,  we  can  mahe  a  first  order 
perturbation  approximation  for  U.  It  is  convenient  to  expand 

+  (*-3) 

k  '  ■'» 

where  the  momentum  k  takes  the  same  set  of  values  as  in  (3.2).  Substituting  this  into  the 
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Schrodinger  equation  (2.3)  then  yields 

ih^a{k,t)  =  f  dt'[coa{Fat' /h  +  k'a)  -  cos{Fat'/h  +  fca))}, 

^  u  ^  •'O 

(5.4) 

where 

V(k  -k')=^'^  ^-iik-k')an  (5  5) 

n 

To  first  order  in  the  potential,  we  have 
a{k,h/{Fa))  =  o(fc.O) 

+  a(k',  0)V(k  -  Jo(sin(()b  -  k')al2)2RI[Fa)), 

xFa^ 

where  Jo(^)  ^  zeroth  order  Bessel  function.  Now,  if  we  make  a  further  expansion, 

a{k,t)  =  J3e-‘*“'”&(m,0,  (5.7) 

m 

where  m  takes  integer  values  1,2,  ..,N  ,  then 

b{m,k/{Fa))  =  6(m,0){l  +  ^  ^  V(fc)e’*“’”Jo(sin(ko/2)2fi/(Fo))}  (5.8) 

Therefore,  U  is  diagonalized  to  first  order  in  V  by  the  transformations  (5.3)  and  (5.7), 
with  the  eigen-frequencies  given  by 

hu;m  =  X;  '^(*)«’*“”‘-^0(Bin(ika/2)2i2/(Fo)).  (5.9) 

k 

This  expression  can  also  be  written  in  terms  of  Vn  as 

fUJm  =  dsJ^^^{zos{27^a)2R/{Fa)),  (5.10) 

where  we  have  taken  the  limit  N  oo. 

When  the  potential  is  zero  (or  constamt),  the  eigenvalues  of  U  are  all  degenerate, 
implying  that  any  state  is  an  eigenstate  of  0.  Therefore,  the  approximation  described  in 
the  last  paragraph  is  really  that  of  a  first  order  degenerate  perturbation.  It  can  be  seen 
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from  (5.9)  or  (5.10)  that  the  degeneracy  is  lifted  in  first  order  by  a  general  potential.  The 
perturbed  eigenstates  can  be  found  by  tracing  back  the  transformations  (5.7)  and  (5.3)  as 


Um(n)  = 

^  k 

=  Y,'^rr^{R/{Fa)), 
I 


(5.11) 


where  the  last  summation  is  over  integers  ensuring  the  boundary  condition  (2.2).  In  the 
limit  of  JV  — ►  oo,  we  have 

«m(«)  —  '^tn-n  [Rl[Fa)),  (5.12) 


a  much  simpler  result.  It  is  noticed  that  the  last  expression  is  also  the  eigenstate  of  the 
Hamiltonian  in  (2.3)  in  the  static  representation  of  the  electric  field,  with  an  uniform  site 
potential  and  an  infinite  system  size. 

The  formula  (5.10)  is  very  convenient  when  the  effect  of  a  few  impurities  is  concerned. 
Consider,  for  instance,  the  case  that  Vn  is  zero  everywhere  except  at  the  origin  n  —  0, 
Then  we  have 

hwm  =  Vq  /  ds72m(co®(2»«’s)2i2/(Fo)).  (5.13) 

In  figure  3a,  exact  numerical  results  for  {hu>n,/Vo}  are  plotted  against  i?/(Fa),  where  we 
have  taken  Vq/jR  =  0.1.  The  above  formula  gives  a  very  good  accotmt  of  the  behavior  of  the 
levels.  In  the  strong  field  limit,  the  levels  are  degenerate  except  for  the  one  corresponding  to 
m  =  0  in  the  above  formula.  As  the  parameter  Rj  (Fa)  is  increased  from  zero,  there  are  two 
levels  peeled  off  from  the  degenerate  center  to  second  order  in  /Z/(Fa),  then  another  two  to 
fourth  order  in  iZ/(Fa),  and  so  on.  This  behavior  of  the  levels  can  be  understood  from  the 
localization  property  of  the  eigenstates  (5.12):  as  the  parameter  i2/(Fa)  is  increased,  the 
localization  length  of  the  states  become  longer,  and  therefore  more  aind  more  sites  feel  the 
existence  of  the  impurity.  At  weaker  field  strengths,  the  levels  which  are  already  outside 
the  degenerate  center  begin  to  oscillate,  accompanying  the  continuous  peeling-off  of  new 
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leveb  from  the  center.  In  figure  3b,  {tiujm/iFa)}  are  plotted  against  R/{Fa),  where  we 
have  taken  again  Vq/R  =  0.1. 

For  a  periodic  potential  of  the  form  (4.21),  a  very  compact  formula  for  the  levels  can 
be  derived  from  (5.9)  as 

hum  =  V  coa{^^m)jQ{sin{irq/p)2R/{Fa)).  (5.14) 

P 

Evidently,  the  levels  btmch  into  p  degenerate  branches,  which  is  a  consequence  of  the 
periodicity  of  the  potential,  and  is  true  even  if  the  potential  is  strong.  Remembering  that 
the  spectrum  defined  in  this  section  corresponds  to  those  modulo  Fa/h  in  the  previotis 
sections,  one  can  read  off  the  degeneracy  from  the  periodic  structure  of  the  spectrum  in 
(5.3).  The  above  formula  has  been  tested  numerically,  and  is  shown  to  give  a  quantitatively 
good  account  of  the  oscillatory  behavior  of  the  levels  in  figures  Id  and  le. 

The  formula  (5.14)  shows  a  further  coalesce  of  the  p  branches  into  one,  for  a  infinite 
discrete  set  of  values  of  the  parameter  R/{Fa)  where  the  Bessel  function  vanishes.  This 
is,  of  course,  generally  true  only  to  first  order  in  the  potential  corresponding  to  the  degree 
of  our  approximation.  It  can  be  shown,  however,  that  for  p  =  2  the  degeneracy  occurs 
infinitely  many  times  as  the  parameter  R/{Fa)  is  varied,  even  if  the  potential  is  strong. 
The  degeneracy  of  the  two  branches  corresponds  to  where  the  two  levels  {hufj/{Fa)}  in 
figures  la  and  Id  have  a  distance  of  unity.  When  all  the  eigenvalues  are  the  same,  the 
evolution  operator  O'  is  proportional  to  the  unity  operator,  and  therefore  every  initial 
state  will  evolve  back  to  itself  after  a  period  of  time,  h/{Fa).  As  has  been  seen  many 
times  in  this  section  and  the  previous  sections,  this  is  the  case  for  the  unperturbed  parent 
Bloch  band.  We  now  see  the  recurrence  of  the  states  at  particular  values  of  the  parameter 
R/{Fa),  if  the  parent  band  is  split  into  two.  The  formula  (5.14)  predicts  an  approximate 
recurrence  of  the  states,  when  the  potential  is  periodic  and  weak. 

6.  Conclusion  and  discussion 
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In  this  paper,  the  effect  of  a  constant  and  uniform  electric  field  on  a  split  Bloch  band 
is  studied  through  a  tight  binding  model.  We  have  shown  how  the  WS  ladders  of  the 
parent  band  is  modified  by  a  splitting  potential,  and  how  the  frequency  spectrum  behave 
as  a  function  of  the  field  strength  sweeping  from  the  weak  to  the  strong  limit.  In  the  case 
that  the  parent  band  is  split  into  a  finite  number  of  sub>bands,  we  have  also  established 
a  new  result  concerning  the  time  evolution  of  the  system.  Further,  we  have  developed  a 
theoretical  framework,  which  exploits  the  symmetries  of  the  system  and  allows  for  easy 
and  clean  numerical  calculations. 

The  general  analysis  has  been  guided  by  the  symmetry  group  of  electric  transla¬ 
tions  [14],  namely  the  gauge  covariant  space-time  translations,  which  exhibit  the  physi¬ 
cal  imiformity  of  the  electric  field.  Three  Abelian  subgroups  of  the  electric  translations 
have  been  tised,  each  emphasizing  certain  aspects  of  the  system.  The  results  obtained  in 
different  ways  are  related  and  complementary. 

The  first  subgroup  emphasizes  the  temporal  uniformity  of  the  system.  The  resultant 
dynamical  symmetry  described  in  section  2  enables  us  to  characterize  the  time  evolution  by 
a  set  of  eigen-frequencies,  and  offers  \is  a  general  and  convenient  calculational  framework  for 
a  system  with  or  without  a  spatial  symmetry.  For  large  systems,  the  eigen-frequencies  are 
the  same  as  the  eigenenergies  (divided  by  h)  of  the  Hamiltonian  in  the  static  representation 
of  the  field.  For  finite  systems,  the  approach  offers  the  convenience  of  eliminating  the 
boundary  effect,  if  bulk  properties  are  the  basic  concern.  Yet,  when  necesssiry,  one  can 
always  insert  an  infinite  potential  somewhere  to  effectively  break  the  ring  geometry  into  a 
simply  connected  one  with  fixed-end  boundary  conditions. 

The  second  subgroup  emphasizes  the  spatial  translational  symmetry  that  the  system 
may  have.  An  immediate  consequence  is  the  periodic  structure  of  the  frequency  spectrum, 
revealing  the  WS  ladders  in  a  natural  way.  An  important  result  has  been  obtained  re¬ 
garding  the  spatial  and  temporal  behavior  of  the  states  in  an  isolated  multi-band  system. 
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There  exists  a  set  of  subspaces  to  which  the  states  return  periodically  in  time;  and  the 
constituent  states  in  each  of  the  subspaces  are  localized  (more  than  exponential)  around 
a  given  unit  cell.  This  result  is  a  generalization  of  the  periodic  recurrence  of  the  states 
derived  from  a  single  band  approximation.  Since  a  localized  state  returns  after  a  finite 
time  interval,  the  state  can  only  travel  a  finite  distance.  In  this  sense  we  say  that  the 
states  are  localized. 

The  following  problem  has  been  studied  in  detail:  how  are  different  set  of  WS  ladders 
are  relatively  plac^  J  in  the  frequency  spectnim,  and  how  do  their  positions  depend  on  the 
field  and  other  parameters?  The  study  is  guided  by  the  analytic  results  obtained  in  the 
weak  and  strong  field  limits  and  in  the  weak  potential  limit.  Exact  numerical  results  for  the 
frequency  levels  have  been  displayed  for  several  cases  of  a  model  periodic  potential.  When 
the  field  is  weak,  the  levels  are  shown  to  behave  according  to  the  prediction  of  the  awiiabatic 
(or  single-band)  approximation.  The  Berry  phases  are  in  general  nonzero  and  shift  the 
levels  to  their  right  positions.  The  interband  couplings  show  up  as  level  repulsions  where 
the  adiabatic  frequency  levels  try  to  cross  one  another.  The  sizes  of  the  gaps  generated  at 
the  intersections  depend  exponentially  on  the  field  and  are  very  small  for  weak  field  values. 
The  gaps  widen  as  the  field  gets  strong,  with  the  spectral  structure  smoothly  transformed 
to  that  in  the  strong  field  limit.  As  expected  from  the  analytical  approximation  appropriate 
for  this  limit,  the  levels  become  those  of  the  parent  band,  individually  shifted  by  the  site 
energies,  showing  the  extreme  localization  of  the  states.  The  initial  slopes  of  the  levels,  as 
the  inverse  field  deviates  from  zero,  are  also  correctly  predicted  by  the  theory. 

The  explicit  time  periodicity  of  the  Hamiltonian  is  exploited  by  the  Floquet  theorem. 
The  appropriate  subgroup  of  electric  translations  is  the  one  with  a  basic  time  step  equal  to 
the  period  of  Bloch  oscillation  of  the  parent  band.  Weak  perturbation  of  the  WS  ladders 
of  the  parent  band  can  be  studied  this  way  and  a  simple  and  general  formula  has  been 
derived  for  the  perturbed  levels.  For  a  weak  periodic  potential,  the  WS  leveb  of  the  parent 
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band  are  shown  to  be  modulated  up  and  down  sinusoidally  as  a  function  of  the  inverse 
field.  The  amplitudes  of  modulation  increase  as  the  field  gets  weak,  while  the  frequency 
of  modulation  approaches  a  constant.  Such  behavior  of  the  levels  is  closely  described  by  a 
very  compact  formula.  Another  interesting  phenomenon  revealed  is  the  approximate  (to 
second  order  in  the  potential)  recurrence  of  an  arbitrary  state,  for  an  infinite  set  of  values 
of  the  field  or  other  parameters.  The  period  of  rec\irrence  is  that  of  Bloch  oscillation  of 
the  parent  band.  For  a  two  band  system,  the  recurrence  is  expected  to  be  exact.  Finally, 
the  effect  of  a  single  impurity  has  also  been  considered  and  analyzed  in  detail. 

Having  summarized  our  basic  results,  we  now  make  a  few  remarks  regarding  the  gen¬ 
erality  and  limitation  of  our  approach  and  the  results.  First,  the  notion  of  electric  transla¬ 
tions  should  be  applicable  to  more  general  systems:  one  dimensional  or  higher  dimensional, 
discrete  or  continuous.  In  this  paper  we  have  limited  our  attention  to  a  system  whose  dy¬ 
namics  is  invariant  under  continuous  time  translations.  What  can  happen  if  the  dynamics 
is  invariant  under  discrete  time  translations  only?  Question  of  commensuration  will  arise 
if  the  dynamics  also  has  a  spatial  periodicity.  In  this  case,  an  Abelian  subgroup  of  electric 
translations  exists  as  a  symmetry  group,  only  when  the  product  of  temporal  and  spatial 
periods  is  a  rational  multiple  of  h/F.  The  situation  u  analogous  to  a  two  dimensional 
electron  system  in  a  periodic  potential  and  a  uniform  magnetic  field. 

Secondly,  most  of  our  results  will  remain  valid  if  our  tight  binding  model  is  modified  to  a 
more  general  form:  with  hopping  amplitudes  nonzero  beyond  the  nearest  neighbors  and/or 
dependent  on  positions,  except  for  the  requirement  of  short  ranged  inter-site  couplings.  For 
example,  the  conclusion  about  the  existence  of  localized  invariauit  subspaces  will  remain 
true,  since  it  was  drawn  from  a  symmetry  consideration  and  from  a  fact  concerning  a  finite 
set  of  first  order  differential  equations.  Another  example  is  the  approximate  recurrence  of 
states  when  the  parent  band  is  weakly  perttirbed. 

Thirdly,  difficulties  will  arise,  however,  for  a  continuous  system  when  all  of  the  inter- 
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band  couplings  are  taken  into  account.  According  to  ref.  4,  the  spectrum  will  become 
continuous  without  any  gaps.  Therefore,  the  validity  of  the  notion  of  WS  ladders  become 
doubtful.  The  problem  does  not  bother  us  too  much,  if  only  a  finite  number  of  bands  (or 
a  finite  spectral  range)  is  concerned,  and  if  these  are  bounded  by  large  energy  gaps.  The 
defense  is  to  invoke  the  usual  wisdom  concerning  the  weakness  of  Zener  tunneling  across 
big  gaps,  and  to  regard  the  results  as  valid  for  long,  but  not  infinitely  long,  time  scales. 
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FIGURE  CAPTIONS 


Figure  1.  Reduced  frequency  levels  {hujj/{Fa),j  =  1,  ..,p}  as  functions  of  R/{Fap)  for 
the  discrete  cosine  potential  (4.21).  (a)  qfp  =  1/2,  V/R  =  0.5;  (b)  qjp  =  1/4,  VjR  =  0.5; 
(c)  q/p  =  2/5,  V/R  =  0.5;  (d)  q/p  =  1/2,  V/R  =  0.1;  (e)  q/p  =  2/5,  V/R  =  0.1.  These 
levels  repeat  periodically  both  upwards  and  downwards  in  the  frequency  spectrum.  As  a 
reference,  the  levels  for  the  parent  band  (V  =  0)  are  horizontal  lines  at  integral  heights. 

Figure  2.  Band  energies  {e/j(s),  /?  =  1,  ..,p}  for  the  adiabatic  eigenvalue  problem  de¬ 
scribed  in  (4.2)  in  the  text.  The  potential  is  Vn  =  V  cos(2s’9/p),  with  V  =  0.5.  The  hopping 
amplitude  is  12  =  1.  The  parameter  q/p  is  (a)  1/2,  (b)  1/4,  and  (c)  2/5  respectively. 

Figure  3.  Frequency  levels  as  functions  of  R/{Fa)  for  a  single  impurity  potential,  with 
Vq/R  =  0.1.  The  levels  have  been  folded  into  an  interval  -0.5  <  huj/{Fa)  <  0.5,  corre¬ 
sponding  to  a  level  spacing  in  the  frequency  spectrum  of  the  parent  band,  (a)  {hufm/Vo} 
are  plotted.  Note  the  Bessel-function  like  behavior  of  the  levels,  (b)  {hu;m/(Fa)}  are 
plotted.  Integral  heights  in  this  plot  correspond  to  the  levels  of  the  parent  band. 
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